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Preface 


The  following  report  is  a  written  version  of  an  abstract  "Coordinates  in  Differential 
Geodesy"  which  was  presented  on  August  23,  1991  at  the  XX  General  Assembly  of  the 
International  Association  of  Geodesy  in  Vienna,  Austria.  The  proceedings  of  the  sessions 
on  "General  Theory  and  Methodology",  Section  FV  were  not  published  by  the  meeting  and 
this  report  represents  a  slightly  amplified  version  of  my  oral  presentation. 

The  abstract  addresses  the  various  roles  of  coordinates  in  theoretical  geodesy: 
determinative  coordinates  in  classical  geodesy,  descriptive  coordinates  in  the 
Marussi-Hotine  formulation  of  differential  geodesy,  and  finally  our  new  notion  of 
permissible  coordinates  in  a  reformulation  of  the  Marussi-Hotine  theory. 
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COORDINATES  IN  DIFFERENTIAL  GEODESY' 


Joseph  Zund 


Summaiy: 

This  report  is  a  written  version  of  an  abstract  presented  in  August  1991,  at  the  XX 
General  Assembly  of  the  International  Associadon  of  Geodesy  in  Vienna.  It  concerns  the 
various  roles  played  by  coordinates  in  differentiai  geodesy,  and  proposes  a  new  approach 
which  offers  the  promise  of  significantly  extending  the  ideas  of  Maiussi  and  Hotine. 
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1.  introductiQn 

Differential  geodesy  is  concerned  with  the  application  of  diffeiendal-geometiic 
techniques  to  the  curves  and  surfaces  which  occur  in  geodesy.  It  was  originally  devised 
by  Antonio  Marussi  and  Martin  Hodne  to  give  a  three-dimensional  formuladon  of 
geodesy,  and  to  clarify  its  conceptual  foundadons,  e.g.  to  avoid  the  ab  initio  dependence 
■—of-  its  fundamental  variables  on  a  pardcular  ellipsoidal  reference  system.  In  effect, 
Marussi  and  Hodne  sought  a  leformuladon  of  classical  geodesy  which  would  re-establish 
it  as  a  branch  of  mathemadcal  physics.  Marussi  laid  down  the  basic  mathematical  and 
physical  requirements  of  the  theory  (the  Marussi  condidons),  he  discovered  the  tensor  of 
gravity  gradients  (the  Marussi  tensor),  and  introduced  his  local  astronomical  coordinate 
system;  while  Hodne  gave  a  comprehensive  presentadon  of  the  theory  in  his  treadse 
(Hodne,  1969).  The  seminal  source  mamrial  which  clearly  indicates  their  ideas  is 
contained  in  the  monographs  (Marussi,  1985)  and  (Hodne,  1991). 

A  basic  problem  in  differendal  geodesy  -  a  problem  which  it  inherits  from 
differendal  geometry  -  is  the  selecdon  of  appropriate  local  coordinate/reference  systems. 
Such  a  selecdon  is  dependent  on  both  mathemadcal  and  physical  consideradons,  and  since 
Hodne  made  an  ambidous  attempt  to  solve  it  in  his  treadse,  it  is  called  the  Hotine 
Problem,  (Zund,  1990).  This  problem  is  unsolved,  and  indeed  it  appears  as  difficult  to 
properly  pose  as  it  is  to  solve  it  in  pracdce.  This  abstract  may  be  regarded  as  a  progress 
report  on  our  attempts  to  understand  and  properly  formulate  the  Hotine  Problem.  This  has 
required  a  re-thinking  of  the  role  of  coordinates  in  differendal  geodesy  and  is  based  on  the 
observadon  that  logically  such  systems  fall  into  two  general  classes  which  we  call 
determinative  and  descriptive  coordinates,  respecdvely.  Hence,  coordinates  belonging  to 
different  classes  do  not  have  the  same  properdes  and  quesdons  posed  for  one  class  need 
not  be  meaningful  for  the  other  class.  In  the  following,  for  brevity  we  denote  by 
(r  s  1,  2,  3)  an  arbitrary  coordinate  system  which  may  belong  to  either  class. 
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2.  Detemunarive  Coordinates 

Determinative  coordinates  are  familiar  from  classical  geodesy  where  they  play  an 
active  role  in  the  determination  of  the  geopotendal  N  as  a  solution  of  the  Modified 
Laplace  equation  in  a  rotating  system: 

AN  =  -2S^  (1) 

wliere  A  is  the  3-dimensional  Laplace  operator  and  m  is  the  (constant)  angular  velocity 
of  the  Earth.  In  (1)  the  unknown  function  N  is  a  dependent  variable,  and  the  chosen  x 
are  independent  variables.  The  choice  of  the  directly  affects  the  prospects  of  solving 
(1)  for  N.  In  effect,  N,  is  as  general,  or  special,  as  the  coordinates  employed  in 
determining  it  An  immediate  concern  in  solving  (1),  or  for  that  matter  any  of  the  partial 
differential  equations  occurring  in  classical  geodesy,  is  whether  the  solutions  are 
susceptible  of  a  meaningful  geometrical/physical  interpretation.  This  is  dramatically 
illustrated  by  the  formal  general  solution  of  the  Laplace  equation 

AN  =  0  (2) 

in  a  Canesian  coordinate  system  y*^  *  (x,  y,  z): 

N  =  f(ax  +  by  +  cz) , 

where  f  is  an  arbitrary  smooth  function  and  a,  b,  c  are  arbitrary  complex  constants 
subject  to  the  condition 

a^  +  b^  +  c^  =  0  . 

This  solution  is  global,  since  y*^  is  a  global  coordinate  system  on  3-dimensional 
Euclidean  space  E^,  however  it  has  no  immediate  geometrical/physical  significance. 

3.  Descriptive  Coordinates 

In  differential  geodesy  as  conceived  by  Marussi  and  Hotine,  the  function  N  is 
assumed  to  be  known,  and  the  analysis  is  then  directed  to  describing  the  geometry  and 
physics  of  the  given  geopotential  field.  In  this  case  the  x*^  play  a  passive  role  whose  only 
purpose  is  a  descriptive  one;  hence  such  x^  may  be  called  descriptive  coordinates.  Then 
(1)  becomes  an  identity,  or  as  one  might  more  appropriately  say  when  concerned  with 
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describing  the  geometry  and  physics  of  the  N-field,  a  consistency  equation  which  couples 
N  with  the  geometric/physical  properties  of  the  field.  This  geometry  is  formulated  in 
terms  of  a  3-dimensional  reference  system,  called  a  vectorial  3-leg,  whose  first  and  second 
vectors  are  associated  with  tangents  to  the  equipotential  surfaces,  while  the  third  vector  is 
normal  to  these  surfaces  and  consequently  tangent  to  their  plumblines.  The  dual  reference 
■— ^stem  is  a  Pfaffian  3-leg  and  hence,  the  entiie  situation  can  be  described  in 
differential-geometric  language  without  the  a  priori  selection  of  a  particular  local 
coordinate  system.  Indeed,  by  use  of  the  leg  calculus,  one  resolves  all  quantifies  along  the 
leg  systems  and  so  the  principal  variables/quanfities  appearing  in  the  theory  are 
coordinate-free.  Since  these  quantifies  depend  only  on  the  choice  of  the  3-leg  system  they 
are  .called  leg  coefficients.  The  term  leg  is  a  literal  translation  of  the  German  word  Bein 
and  was  introduced  by  the  Viennese  geometers  Duschek  and  Mayer  (1930).  It  has  been 
frequently  employed  in  mathemafical/theorefical  physics,  e.g.  in  general  relativity,  and 
more  recently  Grafarend  (1986)  has  suggested  its  use  in  differential  geodesy. 

For  example,  it  can  be  shown  that  (1)  is  replaced  by  the  Bruns  equation 

"ot  ■  ™/3  ® 

where  "/"  denotes  a  leg  (directional)  derivative  along  the  plumbline,  and  H  is  the  mean 
curvature  of  the  equipotential  surfaces.  Upon  introducing  the  local  gravity  n  :=  (in 
Hotine  notation),  the  above  equation  becomes 

nf^-2Hn=  -20?  (4) 

which  specifies  the  change  in  gravity  along  the  plumblines.  Equation  (4)  is  a  consistency 
requirement  which  the  components  of  the  third  leg  vector  v^,  H,  and  a  must  satisfy.  If, 

3 

as  Hotine  does  throughout  his  treatise,  one  imposes  the  Marussi  Ansatz  x  :=  N  with 
1  2 

(x  ,  X  )  being  local  surface  parameters  on  the  equipotential  surfaces,  then  (4)  reduces  to 
merely 

-  2Hv^  =  -2a?  .  (5) 
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Thus,  the  distinction  between  dependent  and  independent  coordinate  variables  is  blurred, 
and  the  components  of  the  leg  variables  X,  ft,  v  become  the  basic  variables.  The  local 
coordinates  are  related  to  the  contravariant  components  of  the  leg  variables  by  the  leg 
differendal  equations 


The  full  set  of  leg  equations,  which  constitute  the  consistency  requirements  describing 
how  the  physical  quantities  N,  n  and  at  are  coupled  to  the  leg  coefficients  defirung  the 
geometry  of  the  N-field,  ate  called  the  Hotine-Marussi  equations  (although  not  all  of  these 
equations  explicitly  appear  in  their  work,  since  neither  of  them  employed  more  than  a 
rudimentary  version  of  the  leg  calculus).  The  unknown  contravariant/covariant 
components  of  the  leg  vectors  may  be  chosen  as  desired  whenever  they  satisfy  the 
Hotine-Marussi  equations. 

4.  Hotine's  Descriptive  Coordinate  Systems 

In  his  treatise,  Hotine  endeavored  to  solve  the  Hotine  Problem,  i.e.  construct  local 

coordinate  systems  subject  to  the  general  requirements  of  the  Marussi  Conditions  (which 

he  did  not  formally  state  or  endorse).  Based  on  the  Marussi  Ansatz  he  proposed  a 

hierarchy  of  five  systems: 

(OJ,  <p,  N)-system  (Chapters  12,  13,  14); 

1  2 

(x  ,  X  ,  N)  normal-system  (Chapter  15); 

triply-orthogonal  system  (Chapter  16); 

(ta,  >^)-system  (Chapter  17); 

synuneuical  (<B,  0,  /^)-systems  (Chapter  18). 

The  diply-orthogonal  system  may  be  discarded  by  virtue  of  the  falsity  of  the  Hotine 

Conjecture  (see  Zund  and  Moore,  1987),  and  it  nims  out  that  the  (co,  0,  /^)-systems  are 

highly  restrictive  and  can  simply  be  regarded  as  specializations  of  the  (or,  (ft,  N)  and 
1  2 

(x  ,  X  ,  N)  normal-systems.  The  (o),  N)-system  is  essentially  the  local  astronomical 

1  2 

system  of  Marussi,  while  the  (x  ,  x  ,  N)  normal-system  is  entirely  due  to  Hotine.  Both 
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of  these  systems  can  be  characterized  in  terms  of  the  contravariam  components  of  the  leg 


vectors  as  follows: 


=  aK  k\  0) .  /  =  0) ; 


/  *  (v^  v^,  h)  in  the  (to,  N)*system  , 

=  (0,  0,  «)  in  the  (x^  x^,  N)  nonnal-systein. 

1  2 

The  (Q),  p,  N)  and  (x  ,  x  ,  N)  ncnrnal-systems  are  quite  different  and  the  former  may  be 
regarded  as  a  specialization  of  the  latter,  but  not  conversely.  It  is  possible  to  develt^  both 
from  a  more  general  notion  which  I  call  the  sistema  finale,  or  F-system.  Using  this  system 
it  can  then  be  shown  that  the  Marussi  Ansatz  involves  no  loss  of  generality  since  under  a 
local  coordinate  transformation  h  x*^,  one  may  choose 

i3_p 

^  “X  fX  “A* 

with  the  subsequent  ad  hoc  identification  of  the  arbitrary  smooth  function  F  with  the 


geopotential  N.  Thus,  the  more  general  F-system  provides  a  basic  framework  which 
includes  both  the  (cu,  p,  N)  and  (x  ,  x  ,  N)  normal-systems  as  limiting  cases.  It  might 
be  noted  that  in  a  sense  for  descriptive  systems,  the  question  of  the  generality  of  the 
systems  is  a  moot  issue:  a  chosen  x*^  lacks  generality  only  to  the  extent  that  they  fail  to 
provide  a  description,  which  is  convenient,  or  otherwise,  of  a  particular 
geometrical/physical  field  property. 

Thus,  the  Marussi-Hotine  approach  to  differential  geodesy  essentially  contains  two 
descriptive  local  coordinate  systems,  which  as  noted  above  are  structurally  different  For 
example,  the  (tu,  p,  N)-system  is  applicable  only  to  a  single  equipotential  surface,  viz  it  is 
'frozen'  on  the  surface.  More  precisely,  if  one  considers  a  family  I  of  such  surfaces  then 
one  can  choose  the  longitude  and  latitude,  i.e.  (to,  P),  as  local  surface  parameters  on  any 
surface  S  of  £ ,  then  one  can  take  the  same  (ox  P)  as  surface  parameters  on  a 
neighboring  equipotential  surface  only  when 

(i)  S  is  a  developable  surface,  or 

(ii)  the  congruence  F  of  plumblines  of  S  are  straight  lines. 
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Both  (i)  and  (ii)  are  severe  restrictions:  (i)  requires  that  the  Gaussian  curvature  K  of  S 

is  zero,  so  typically  S  would  be  a  cylinder,  cone,  etc.  in  contrast  to  our  expectadons  of 

K  >  0  for  a  sphere,  spheroid  or  ellipsoid;  (ii)  requires  that  the  curvature  %  ^he 

plumblines  vanish  so  T  is  a  linear  congruence,  hence  S  could  be  a  sphere  but  not  a 

spheroid  or  an  ellipsoid!  Hodne  was  aware  of  (ii),  and  devised  his  isozenithal 

differendadon  (see  Chapter  14,  of  (Hodne,  1969))  to  work  around  it  However,  in  view  of 

(i)  such  a  manuever  is  of  dubious  value.  The  frozen  character  of  the  (co,  0,  N>system  is 

an  especially  heavy  impediment  in  Hodne's  work  since  many  of  his  calculations  in  the 

1  2 

treatise  are  given  in  this  system,  or  one  step  away  from  it.  The  (x  ,  x  ,  N) 

normal-system  was  explicidy  constructed  to  provide  a  local  coordinate  system  for  which 
1  2 

(x  ,  X  )  can  be  taken  as  local  parameters  on  the  surfaces  of  I  in  such  a  manner  that 
1  2 

(x  ,  X  )  do  not  change  along  T.  The  notion  of  such  a  system  is  intuitively  appealing  and 

seems  to  have  originated  with  Hodne,  but  is  not  a  standard  piece  of  Gaussian  differential 

geometry.  The  centerpieces  of  his  theory  arc  his  varictional  equations  which  prescribe 

how  the  intrinsic  and  extrinsic  properties  of  the  surfaces  of  Z  vary  along  F.  These 

equations  are  undeniably  elegant  and  a  real  tribute  to  Hotine's  ingenuity.  However,  it  is 

not  clear  whether  they  lead  to  useful  results.  For  example,  they  give  no  clue  as  to  what 
1  2 

the  (x  ,  X  )  are,  or  how  to  choose  them,  or  even  whether  a  non-trivial  pair  of  such 

parameters  exist.  Indeed,  the  variational  equations  are  very  complicated  and  appear  to  be 

integrable  in  closed  form  only  when  £  reduces  to  a  family  of  concentric  spheres,  i.e.  F 

is  a  linear  congruence.  This  need  not  be  a  serious  defect,  since  the  proper  role  of  these 

equations  is  that  of  consistency,  not  determining,  equations;  however,  by  merely 

appending  them  to  the  Hotine-Marussi  equations  one  greatly  increases  the  complexity  and 

1  2 

the  number  of  these  equations.  Thus,  neither  the  {(O,  0,  N)  nor  the  (x  ,  x  ,  N) 

normal-system  is  truly  satisfactory,  or  really  workable.  Both  are  heavily  restricted,  and  it 

1  2 

IS  possible  that  the  more  general  (x  ,  x  ,  N)-system  as  derived  i"rom  our  F-system  theory 
will  prove  more  successful.  What  is  clear  is  that  —  at  present  —  we  do  not  possess  a  really 
non-trivial  solution  of  the  Hodne  Problem. 
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The  purely  detenninative/descriptive  role  of  the  is  extreme:  in  classical 
geodesy  the  geometric  content  of  the  theory  may  be  difficult  to  determine  and  often 
appears  to  be  contrived,  while  in  differential  geodesy  —  in  the  Manissi-Hotine  formulation 
-  one  is  faced  with  devising  a  natural  description  of  the  properties  of  the  geopotendal 
field.  Both  approaches  may  be  termed  'extractive'  theories  where  one  is,  in  effect, 
-  -challenged  to  extract  the  geometry  from  a  physical  situation.  These  formulations  are 
perfectly  valid  and  workable,  but  neither  is  completely  satisfying.  Perhaps  this  is  what 
Hotine  sensed  but  did  not  put  his  finger  on,  when  in  (Hotine,  1965)  he  wrote: 

'Tor  some  time  I  have  had  an  uncomfortable  feeling  that  there  is 
something  missing,  either  in  the  classical  theory  of  gravitation  itself,  or  in 
the  mathematical  handling  of  it” 


Quite  distinct  from  the  Marussi-Hodne  desciipdve  approach  one  may  consider  the 
situadon  in  differendal  geodesy  from  an  alternate  more  robust  viewpoint.  We  call  this  a 
permissible  viewpoint  in  which  one  specifies  certain  of  the  dynamical  and  geometric 
quanddes,  e.g.  the  leg  coefficients:  ^2^  /j.  Jy  *^1’  ®2’  Hodne  (1969)  and 

Zund  (1990),  and  then  asks  whether  the  Hodne-Marussi  equadons  are  compadble  with  this 
choice.  In  this  case,  one  now  seeks  the  geopotendal  fiincdon  N  subject  to  certain 
dynamical/geometiic  condidons  and  (3),  or  (4),  is  no  longer  an  idendty  but  a 
determinadve  equadon  and  the  solution  of  this  equadon  must  be  consistent  with  the  full 
set  of  Hodne-Marussi  equadons  which  include  the  Lamd  equadons  for  the  flamess  of  E^. 
Depending  on  the  prescribed  condidons  some  of  these  Lamd  equadons  may  be  idenddes 
while  others  are  non-trivial  determining  equadons  for  the  unspecified  leg  coefficients.  For 
example,  one  might  seek  a  soludon  for  N  in  which  K  :=  <^2  *  ^  ® 

one  having  equipotendal  surfaces  of  non-negadve  Gaussian 
curvature  and  curved  plumblines,  subject  to  certain  non-zero  components  of  the  3-leg 
vectors  A,  fi,  and  v.  Given  such  a  specificadon  either  the  Hodne  Problem  is  solvable  or 
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insoivable,  and  both  situations  are  interesting.  Relative  to  the  former,  one  has  constructed 
a  geopotendal  field  having  certain  prescribed  dynamical/geometric  properties;  while  in  the 
latter  case,  assuming  that  the  resulting  Hodne-Manissi  system  of  equadmis  are  tractable, 
one  may  show  that  the  prescribed  conditions  are  inconsistent  or  that  the  conditions  are 
inadequate,  i.e.  not  sufficiendy  specified,  to  permit  a  soludon.  Obviously  the  former  case 
is  the  most  interesdng,  but  the  latter  is  valuable  in  that  it  may  indicate  how  the  Hodne 
Problem  should  (or  should  not)  be  formulated.  This  approach  is  said  to  be  permissible, 
since  it  involves  exhibiting  only  those  soludons  of  N  and  the  unspecified  leg  coefficients 
which  are  permitted  by  the  assumed  dynamical/geometric  condidons.  The  final  step  of 
such  a  procedure  would  be  to  solve  (6)  for  the  resulting  permissible  coordinates  tJ.  Such 
a  methodology  would  be  an  active  rather  than  a  passive,  or  'extracdve',  approach  to 
differendal  geodesy,  and  could  lead  to  a  more  robust  and  vigorous  formuladon  of  the 
geometric  aspects  of  theoredcal  geodesy. 

Thus,  we  believe  that  differendal  geodesy  is  a  live  and  vibrant  theory,  a  theory 
which  may  well  step  beyond  the  goals  set  for  it  by  the  visionary  efforts  of  Marussi  and 
Hodne.  Reladve  to  it  and  its  prospects  for  the  future  I  would  concur  with  the  famous 
comment  of  Sir  Winston  Churchill: 

"Now  this  is  not  the  end.  It  is  not  even  the  beginning  of  the  end.  But  it 

is,  perhaps,  the  end  of  the  beginning." 


8 


REFERENCES 


Duschek,  A.  and  Mayer,  W.  (1930).  Lehrbuch  der  Diffemdalgeometrie  Band  II 
Riemannsche  Geometric,  B.G.  Teubner,  Leipzig. 

Grafarend,  E.  (1986).  Differential  geometry  of  the  gravity  field,  Mansucripta  Geodaetica 
11,29-37. 

Hodne,  M.  (1965).  Trends  in  mathemadcal  geodesy,  Bolledno  di  Geodesia  e  Scienze 
Affini,  anno  JOOV,  607-622;  reprinted  in  Hodne  (1991)  as  Paper  1,  5-21. 

Hodne,  M.  (1969).  Mathemadcal  geodesy,  U.S.  Department  of  Commerce,  Washington, 
D.C. 

Hodne,  M.  (1991).  Differendal  geodesy,  edited  with  commentaries  by  J.  Zund, 

Springer- Verlag,  Berlin. 

Marussi,  A.  (1985).  Intrinsic  geodesy,  translated  by  W.L  Reilly,  Springer-Verlag,  Beriin. 

Zund,  J.  (1990).  The  Hodne  problem  in  differendal  geodesy,  Manuscripta  Geodaedca  15, 
373-382. 

Zund,  J.  and  Moore,  W.  (1987).  Hodne's  conjecture  in  differendal  geodesy.  Bulletin 
G6od6sique  61,  209-222. 


9 


